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ABSTRACT 

Context. In this paper we study the dynamics of the v Andromedae planetary system proposed by Curiel et al. (201 1). We focus on the 
study of the 3:1 Mean Motion Resonance (hereafter MMR) between u Andromedae-d and the recently discovered v Andromedae-e 
(hereafter v And-d and u And-e). 

Aims. Numerical simulations of the dynamics of the four planet system are conducted. 

Methods. The previously reported apsidal resonance between v And-c and v And-d is confirmed. In addition, we find that u And-d 
and v And-e are also in an apsidal resonance condition. 

Results. Our results further indicate that the u Andromedae planetary system configuration is in the middle of a stability island in the 
semimajor axis-eccentricity domain. Additionally, we performed numerical integrations of the planetary configuration over 500 Myr 
and found it to be stable. 

Conclusions. We conclude that, within the uncertainties in the value of the orbital parameters, it is likely that v Andromedae planetary 
system will remain stable for a long timescale. 

Key words, planets and satellites: dynamical evolution and stability -stars: individual: v And 



1. Introduction 

About ~ 10% of the exoplanets known to date are in multiple 
planet systems (http://exoplanet.eu). This is likely to be a lower 
limit since there is still a strong observational bias towards the 
detection of planets with short periods. Therefore, some of the 
extrasolar planets detected so far may ha ve planetary compa n- 
ions not yet detected on distant orbits (e.g. lCorreia et al.l d2009l) ). 

Of the exoplanets discovered in multiple planet systems, sev- 
eral cases have been found to be in mean motion resonances 
(MMR). This is important beca use resonances tend to stab ilize 
the orbits of the planets involved. lMurrav & Dermotg(l2000t) sec- 
tion 8.3. Of the cases discovered to date, the 2:1 mean motion 
resonance is the most common (e.g. HD 73526, HD 82943, HD 
128311, GJ 876, Kepler 9 and HD 37124), but there are other 
configurations such as the the 3:2 MMR in HD 45364, the 3:1 
MMR in HD 75732 and the 5:1 in HD 202206. In our Solar 
System, the 3:2 MMR between Neptune and Pluto (albeit not 
a bona-fide planet) that allows the latter to be in a high ec- 
centricity and high inclination on an orbit that crosses the or- 
bit of Neptune, su ch that these two objects will never collide 
(e.g. lVaradH d!999l) ). In the context of the so-called Nice model, 
it is believed that in the early Solar System, Jupiter and Saturn 
crossed the 2: 1 resonance, leading a significant re arrangement of 
the ge n eral architecture of our plan etary system ( iTsiganis et al.1 
d2005h . lMorbidem et al.l (12001 and lGomes et "all (120051) ). 

It has been pointed out recently bv lCorreia et al.1 d2009t) . that 
the study of the dynamics of mean motion resonances of two 
or more planets interacting in a system, offers the opportunity 
to constrain and understand the process of planetary formation 
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and evolution, since these resonances most probably arise from 
planetary migration. 

The v Andromedae (hereafter v And) planetary system was 
the first multiple, ext rasolar planetary s ystem discovered orbit- 
ing a solar-type star dButler et al.l 1 1999b . The system is known 
to harbour three extrasolar planets with masses ranging from 
0.69 to 14.57 Mj, with Mj being the mass of Jupiter. Recently, 
ICuriel et"ai] d201 lb have found a fourth planet orbiting v And 
named, as is the convention, u And-e. It is the pourpose of the 
work presented in this paper to analyse the dynamics of this new 
planet to determine whether it is in MMR with v And-d and, 
hence, in a likely stable configuration. 

The paper is organized as follows. In Section 2 we review the 
general characteristics of the central star v And and the known 
orbiting planetary system. In Section 3, we perform a global 
frequency analysis of the two planets possibly involved in the 
3:1 resonance, namely, v And-d and v And-e. Additionally, we 
present results on the long term evolution of the system. Section 
4 focuses on the apsidal resonance between v And-c and v And- 
d and Section 5 focuses on the apsidal resonance between v 
And-d and v And-e. Finally, we summarize our results and offer 
our conclusions in section 5. 



2. v Andromedae system characteristics 

v And is a bright F8 V star with a mass of 1.3 M and stel- 
lar radius of 1.56 R Q dButler et al.lll999h . The distance to the 
star is estimated to be about 13.47 pc, that is 43.93 lyrs 
( Perrvman et alj T997I) . The estimated age of the star is 5 Gyrs 
dBaliunas et al]| 19971). and its r o tational per i od is b etween 9 and 
12 days dBaliunas et al.l dl997l) . iFord et al.1 dl999h ). v And was 
the first multiple exoplanetary system detected around a main 
sequence star. It was reported as being a triple planetary sys- 
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Tabl e 1. Properties of the planetary companions in v And (from 



lable 1. Properties < 
ICuriel et all (hoi 11) 1 



Planet 


M(Mj) 


a (AU) 


e 


b 


0.6876(44) 


0.05922166(20) 


0.02150(70) 


c 


1.981(19) 


0.827774(15) 


0.2596(79) 


d 


4.132(29) 


2.51329(75) 


0.2987(72) 


e 


1.059(28) 


5.24558(67) 


0.00536(440) 



Table 2. Fundamental frequencies calculated from the nominal 
solution 





Frequency (°yr) 


Period (yr) 


n d 


0.285194 


1262.3 


n e 


0.0947104 


3801.06 


ft 


0.0000257259 


38871.3 


ge 


0.00109351 


914.484 


<% 


1.52241 


236.467 




0.415535 


866.353 



tem by iButler et"aT] (Il999h . The estimated masses and orbital 
paramenters for each of the planets around v And are listed in 
Table Q] Several studies have been made about the origin and 
dynam i cal stability of the triple planetary system ( Ford et al.l 
d2005h . lCMang et al.ld200ll) . IRivera & Haghighipourl d2007l) and 
references therein). The main conclusion of these st udies, con- 
siderin g only the first three planets discovered are: iFord et all 
(120051) found out that planet-planet scattering with a fourth lost 
planet co uld explain the high eccentricities of v And-c and 
v And-d. IChiang et alJ d200ll) found out that the apsidal res- 
onance between u And-c and v And-d is observed for mu- 
tual inclinations, betwe e n thes e two planets, smaller than 20°. 
iRivera & Haghighipourl (120071) concluded that there are stabil- 
ity regions for test particles around the 3:1 and 5:1 MMR, and 
that for a > 7.5 AU all test particles were stable for at least 
10 7 yrs. They reported additionally that test particles just out- 
side the 1:3 MMR with v And-d experience large oscillations 
reaching eccentricities up to the range of 0.2 to 0.3; but these 
particles are protected from cl ose approaches with v And-d by 
the e-w mechanism. Recently. ICuriel et al.1 (1201 ll) have discov- 
ered a fourth planet orbiting the system on the basis of a refined 
fit for the radial velocity data . Its properties are also described 
in Table fTI ICuriel et aD d201 ll) propose that the system is close 
to a 3:1 MMR resonance. The initial estimate for the period of 
the fourth planet (u And-e) is 3848.86 days, and the period of 
the third planet {u And-d) is 1276.46. The ratio between the two 
periods is 3. 02, very close to being in exact resonance. In a pre- 
vious study, IRivera & Haghighipourl (120071) find that there is an 
island of stability in the semimajor axis-eccentricity parameter 
domai n, just outside the externa l value s corresponding to the 3:1 
MMR. IRivera & Haghighipourl (120071) use a large collection of 
test particles to sample the possible location of new planets. 

3. Orbital stability 

In this section we analyse the dynamical stability of a planetary 
system characterized by the orbital parameters reported in Table 
Q] We refer to this set of parameters for the v And planetary sys- 
tem as the nominal solution. There are different procedures to 
check the stability of a sy stem using the frequency analysis de- 
veloped by |Easkai] ^199^J^^^^i^^^M^}^sfoUow 
the approach taken bvlCorreia et al.l (120051) . ICorreia et alJ d2009l) 
and ICouetdic et al.l (1201 Oh to analyse the stability of the plane- 
tary configuration. 
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Fig. 1. Time evolution of the three resonant angles for the 3:1 
resonance between v And-d and v And-e. The angles 9\ and 8o 
circulate, while #3 librates. The bottom panel shows the oscilla- 
tion in the time evolution of the eccentricity of u And-e (called 
here ei) 



3.1. The 3:1 resonance 

To test whether the v And system is trapped or not in the 3:1 
MMR resonance, we perform a frequency analysis of the nomi- 
nal solution computed over 10 Kyrs. Since the orbit of the inner 
planets of the system is well constrained, we fix the orbital pa- 
rameters of v And-b, v And-c and v And-d according to Table 
Q] We also assume the orbits of all planets are coplanar. A series 
of initial conditions for the eccentricity (e) and semi-major (a) 
axis of planet e are constructed to test the stability of the region 
around the nominal solution. 

For each initial condition, the orbits of the planets are in- 
tegrated over 10 kyrs using the hybrid integrator included in 
the Mercury 6 code dChambersI 1 1 9991) . For most of the inte- 
gration, Mercury uses a mixed-variable symplectic integrator 
(IWisdom & Holmanll99ll) with a time step approximately equal 
to a fiftieth (wl/50) of the Keplerian orbital period of the closest 
planet (in this case u And-b). During close encounters, Mercury 
uses a Bulirsch-Stoer integrator with an accuracy parameter of 
10~ 12 . The stability of the orbit corresponding to each set of (a, e) 
values i s mea sured using the freq uency analysis introduced by 
lLaskarl (1 19901) and lLaskarl dl993l) . According to lLaskarl (119931) 
the difference (D) in the value of the fundamental frequency of 
the motion of the planet under consideration, obtained over two 
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Fig. 2. Stability analysis of the v And system as a function of semi-major axis, a, and eccentricity for planet-e, e. The logarithm of 
the stability parameter D is shown on a gray scale map. The lighter colors correspond to the more unstable systems. The log-value 
+2 for the parame ter D corresponds to the ejection of planet-e. The ir symbol corresponds to the nominal solution reported by 
ICuriel et"a"D d201 ll) . Panel b) is a close-up, in a, of the wider domain shown in panel a) of the parameter space around the nominal 
solution. 
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Fig. 3. Orbits in rotating reference frames for 500 yrs for the v And system. Panel (a) shows the reference frame with the angular 
velocity M c in which the 3:1 MMR resonance is observed. Panel (b) shows the same points but in a reference frame that is rotating 
with Md, where we observe the 1:3 MMR resonance 



consecutive time intervals is a measure of the secular stability of 
the trajectory. 

In addtion, we also identify strongly unstable systems as 
those in which: 1) two planets collide, 2) a planet hits the star (if 
its astrocentric distance is > 0.005 AU), or 3) a planet is ejected 
from the system (assumed to occur if the planet travels beyond 
an astrocentric distance of 100 AU). 

The thre e possible resonant argume nts of the 3:1 resonance 
according to lMurrav & Dermottl (l2000l) p.491. with Q = 0.0 be- 
cause we consider coplanar orbits, are the following: 

6 1 =3A e -A d -2co d , (1) 
6*2 =3Ae- A d - to e - a) d (2) 



and 

#3 = 3A e - Aj - 2at e , (3) 

where Ad, a>d and A e , a> e are the mean longitude and argument 
of the pericenter of planet e and d, respectively. Fig. Q] shows 
the time evolution of the three possible resonant angles and the 
time evolution of the eccentricity of planet-e, in order to make 
a direct comparison of its periodicity with the resonant angles. 
We find that 6\ and 02 are circulating while 63 is librating around 
zero. This behaviour means that the two planets are closer to the 
resonance defined by #3. 

It is clear form Figure[TJthat 02 and 63 have periodic and sud- 
den changes with periods of 866.353 yrs. (we found this period 
using a Fourier analysis), this period coincides with the period of 
approximately 866.353 years that can be seen in the oscillation 
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Fig. 4. Apsidal resonance between planet-c and planet-d. The 
two eccentricities are anti-correlated as shown in the top panel. 
The eccentricity of planet c (here ei,) gets close to zero periodi- 
cally. We show in the middle panel the argument of the pericen- 
tre of the two planets, it can be noticed that they tend to follow 
each other. Finally, in the bottom panel we show the difference 
between the arguments of pericentre (Aa>23) of planets c and d. 
As we can notice this angle is librating around zero. 
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Fig. 5. Apsidal resonance between planet-d and planet-e. The 
two eccentricities of planets d and e (shown here in light grey and 
dark grey, respectively) are anti-correlated. The middle panel 
shows the arguments of pericentre of these two planets, it is pos- 
sible to notice that they follow each other. Finally, in the bottom 
panel is shown the difference between these two arguments of 
pericentre (AW34), it is possible to notice that this angle is librat- 
ing around zero. 



of the eccentricity of planet-e about zero. This behaviour can be 
explained by the fact that to e is not well defined for circular or- 
bits. Therefore #3 is not circulating every 866.353 yrs., but rather 
librating and every period that planet e comes back to zero this 
resonant angle suffers the sudden changes observed in Fig.Q] 

The fundamental frequencies of the system are the two mean 
motions (known also as mean angular velocities) rid — ^jf and 
n e = the two secular frequencies of the pericenter gj = 
^f- and g e = and the libration frequency of the resonant 

argument lg — -Jr. The values of each one of these are shown in 
Table 

3.2. Stability analysis 

To analyse the stability of o ur nominal so lution, we perform 
a global frequency analysis dLaskarl [l993) in its vicinty (Fig. 
[2]> (a) and (b), similarly as it h as been done bv ICorreia et alJ 
d2005l) . ICorreia et all (120091) and lCouetdic et al.1 'fcOlOh . For the 
two planets that are possibly involved in the the resonance (d 
& e), the system is integrated using a two-dimensional mesh of 



11,165 initial conditions. We vary the semi-major axis and ec- 
centricity of planet e, and keep the orbital parameters of the other 
planets at their nominal values. 

We integrate each case numerically over 10 kyrs, and then 
we compute the stability indicator D. This stability indicator is 
computed by comparing the variation in the measured mean mo- 
tion over two consecutive T — 5 kyr intervals. The parameter D 
is defined as: 

\n e — n'\ 

D = (4) 

where n e and n' e are the mean motion of planet e in the first 5 
kyr and last 5 kyr of the integration, respectively; the units of D 
are deg/yr 2 . This parameter D is a measurement of the chaotic 
diffusion of the trajectory. It should be close to zero for a reg- 
ular solution and it has a high value for strong chaotic motion 
(lLaskarll993l) . We also performed 10 integrations of 1 Myrs and 
found that in the v And system, a value of D < 10~ 7 is required 
for regular motion (comparing the D of stable and unstable inte- 
grations). 
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The left panel of Figure[2](a) shows the results of our mesh of 
11,165 numerical simulations in the domain of semimajor axis, 
a 6 [4.7,5.8] AU, and eccentricity e e [0,0.5]. From Fig.|2](a) 
we see that there is an island of stability around the nominal so- 
lution, shown here with the symbol ~k (located at a = 5.24558 
AU). The region does not seem to contain ejected particles (des- 
ignated with D — +2 in this plot). We show also the location 
of the 3:1 MMR as a vertical line at a = 5.22672 AU. At this 
scale the location of the nominal solution and the location of the 
3:1 MMR are very close to each other. Therefore we decided to 
zoom-in around the region of the nominal solution, and explore 
it numerically. 

We performed another set of 11,615 integrations around 
the stability island, focusing on the parameter domain a e 
[5.18,5.30] AU and e e [0,0.5]. Our results are shown in the 
right panel of Fig.|2](b). The shape of the stability island around 
the 3:1 resonance is easily appreciated on this scale. We notice 
that most of the systems explored in the a-e mesh are unstable 
above e = 0.1. The nominal solution reported in ICuriel et al.l 
( 120111) has a = 5.24558 AU and e = 0.00536 (shown again 
with That puts the system right in the middle of the stabil- 
ity island. This resu lt support the idea that the solution found by 
ICuriel et al l d20TTb is stable and robust, as suggested by these 
authors. Also, it appears that the only stable zone that exists in 
the vicinity of the nominal solution is the zone protected by 3:1 
resonance. 

In Fig. [3] (a and b) we show the evolution of v And-d and 
v And-e over 500 yr in the rotating reference frame of the in- 
ner and outer planet, respectively. Since the system is so close to 
the 3:1 resonance the relative position of the two planets is re- 
peated, and their minimum distance in our simulation was 2.47 
AU. The paths of the two planets in this rotating reference frame 
shows the relationsh ip between the resonance a nd the frequency 
of conjunctions (see iMurrav & Dermotj d2000l) . p. 325). In this 
particular reference frame, every three orbits of v And-d corre- 
sponds to one orbit of v And e. Figure[3]also shows the libration 
of each planet around its equilibrium position. Since there are 
two other planets on the system v And, the equilibrium position 
precesses and does not maintain its path, in the rotating reference 
frame, for a much longer time. 

3.3. Orbital evolution 

The previous stability analysis supports the idea that the u And 
system is trapped in the island-like stability zone associated with 
the 3:1 resonance, and that it is very likely to be stable on long 
term evolution. To test this suggestion we carried out a long 
term numerical integration of the nominal solution. Using the 
Mercury 6 code with the same parameters as described above, 
we integrated the orbits of the four planet s for 500 Myrs. W e 
found that, as suggested by the study of ICuriel et"aTI (1201 ll) . 
these orbits are stable in the long term. 

Due to the intense resonant gravitational interaction the three 
planets involved (c, d & e) exhibit significant changes in their 
eccentricities, as can be seen in Figs. |4] and [5] 

4. Apsidal resonance between v And-c and v And-d 

In this section we report on the apsidal resonan ce that was 



observ ed 
(1201 ll). 



in previous re search o n v And (e.g. 



Lissauer & Rivera 
Barnes & Ouinnl 
Nagasawa & Linl (120051) . etc). 



(12001 


). IChians & Murravl 


(12004. 


Michtchenko & Malhotra 



Curiel et al.l 
<|2000h. 
(|2002h . 

(2004), 



The secular apsidal resonance is shown in Fig. [4] We find 
that the eccentricities e d and e\ are anti-correlated as well, that is, 
when e d is maximum e± is minimum and vice versa. We show cj c 
(dark gray) and u) d (light gray) in the middle panel of Fig. [4] The 
difference between these two angles is shown in the bottom panel 
of Fig.|4] it shows that Aco cd oscillates around zero with an initial 
amplitude of ~ 77° and a short period of 6 x 10 3 yrs. Therefore 
showing that planets c and d are indeed in apsidal resonance. 

The physics of the apsidal resonance is the following. The 
eccentricity of v And-c and and v An d-d are quite large (e c = 
0.2596 and e d = 0.2987, Table Q] of (ICuriel et al.l l201lb ) and 
their semi-major axis are not too different from each other (a c = 
0.827774 AU and a d = 2.51329 AU). In principle, they could 
be as close to each other as D m \ n = a d (l - e d ) - a c (l + e c ) = 
0.719906 AU, but this never happens. As it is possible to notice 
in Fig. [4] when the two orbits are aligned (that is Aa) c d = 0.0) 
the eccentricity of planet-d is maximum, but at the same time the 
eccentricity of planet-c is reduced to values e c ~ 0. This helps to 
avoid close encounters between the two planets, hence helping 
to stabilize the system, then the two planets never get closer than 
0.9 AU. 



5. Apsidal resonance between v And-d and v And-e 

In addition to planet-c and planet-d being in apsidal resonance, 
we find that planet-d and pl anet-e are in apsida l resonance as 
well. This was suggested by ICuriel et"aTI (1201 ll) . and here we 
give the details about this resonance for the first time. Figure [5] 
shows this apsidal resonance in action. We find that e d and e\ are 
anti-correlated. We also find that the maximum eccentricity that 
planet-e can achieve is quite big (e max = 0.233105). Therefore, 
in principle, it is possible that planet-e and planet-d can get 
close to each other D m ; n = 04(1 - e\) - a d {\ + e d ) = 2.051726 
AU, but owing to the apsidal resonance, they never get closer 
thanD min = 2.47421 AU. 



6. Summary and conclusions 

We have carried out a stability analysis for the 4 -planet system 
around v And by means of a frequency analysis lLaskarl (1 19931) 
and long term numerical simulations. On the basis of our results, 
we find evidence that v And-d and v And-e are in a 3:1 mean 
motion resonance, as suggested by prev ious studies. The nomi- 
nal solution found bv lCuriel et al.l(120Tll) is located in the middle 
of a island of stability in the a-e paramenter space. 

The v And system is found to be a rich dynamical system, 
with three of the planets discovered to date interacting strongly 
via resonances. As described in the paper, v And e interacts via 
the 3:1 MMR with v And d, and additionally they are in apsidal 
resonance. This prevents them from having close encounters. On 
another hand, the apsidal resonance reported by previous authors 
between v And-c and v And-d is still present when u And-e is 
considered. This means that the three planets are interacting via 
apsidal and MMR resonances; 

The nominal solution of ICuriel et al.l (1201 lh is here proved 
to be stable for 500 Myrs, this along with our results of the sta- 
bility analysis done via the global frequency analysis, allows us 
to conclude that the nominal solution is robust and stable. 

Finally the reported eccentricity of u And-e is at the moment 
very close to zero (<? = 0.00536). However, according to our re- 
sults, it should increase to e — 0.016 in approximately 10 years. 
If it is possible to measure this variation after such period, it will 
provide an important confirmation of the dynamical properties 
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described in this work. Our results further predict that the eccen- 
tricity of v And-e will reach its maximum value of 0.2 in around 
450 years. 
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